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Specific heat per particle (c„) of an ideal gas, in many occasions, is interpreted as energy fluctua¬ 
tion per particle (Ae 2 ) of the ideal gas through the relation: Ae 2 = kT 2 c v , where k is the Boltzmann 
constant and T is the temperature. This relationship is true only in the classical limit, and deviates 
significantly in the quantum degenerate regime. We have analytically explored quantum to classical 
crossover of this relationship, in particular, for 3-D free Bose and Fermi gases. We also have explored 
the same for harmonically trapped cases. We have obtained a hump of Ae 2 /fcT 2 c[, cl ' ) around the 
condensation point for 3-D harmonically trapped Bose gas. We have discussed the possibility of 
occurring phase transition with discontinuity of heat capacity from existence of such a hump for 
other Bose and Fermi systems. 

PACS numbers: 03.75.Hh, 03.75.Ss, 67.85.-d 


I. INTRODUCTION 


Ultracold atomic gas in an optical trap is a favourite 
hunting ground for theoreticians and experimentalists [1- 
3]. Within the last two decades numerous works were 
done on this topic [1-3]. Quite a few major branches 
e.g. Bose-Einstein condensation in harmonic trap [1], 
ultracold Fermi gas in harmonic trap [3], ultracold gas 
in optical lattice [2] , superfluid-Mott insulator transition 
in optical lattice [2, 4], etc have evolved from this sin¬ 
gle topic showing its potential. However, the thermo¬ 
dynamic properties of the ultracold atomic gases have 
always been the centre of attraction. Bose-Einstein con¬ 
densation fraction [5-8], temperature dependence of en¬ 
ergy and specific heat of ultracold Fermi gases [9-11], mo¬ 
mentum distribution for harmonically trapped Bose gas 
[12], momentum distribution for harmonically trapped 
Fermi gas [13], temperature dependence of the chem¬ 
ical potential [14], temperature dependence of critical 
number of particles for collapse of attractively interact¬ 
ing Bose gas [15], temperature dependence of thermody¬ 
namic properties of a unitary Fermi gas [16], etc have 
already been studied in this regard. As the discontinuity 
of specific heat (at constant volume and average parti¬ 
cle number) confirms about occurrence of the second or¬ 
der phase transition, measurement of the specific heat is 
very important. Surprisingly, specific heat for harmoni¬ 
cally trapped Bose gas, as far as we know, has not yet 
been precisely measured, specially around the condensa¬ 
tion point. This is because of difficulties associated with 
the inhomogeneity of the trapped condensate, and the 
corresponding difficulties in analyzing the density pro¬ 
files within local density approximation [17]. If this is so 
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difficult to be measured, is there any other measurable 
thermodynamic variable which confirms discontinuity of 
specific heat? This paper finds a possible answer to this 
question from a general perspective. 

Specific heat is very often interpreted as energy fluc¬ 
tuation. Energy fluctuation, for a single particle subsys¬ 
tem, is nothing but the variance of energy Ae 2 = e 2 — e 2 , 
where e = J2i e iPP e 2 = J2 t e iPii C is the energy eigen¬ 
value of the particle at the ith eigenstate, and pi is the 
thermal equilibrium probability of the ith eigenstate. As 
long as we take this probability to be Pi = e~ ei ^ kT /Z , 
where Z = JA e _e< / fcT j s the canonical partition func¬ 
tion, we obtain the familiar relationship: Ae 2 = kT 2 c£ V> 
where c£ c1 ^ is the classical specific heat of a single parti¬ 
cle. The same probability is eventually applicable for a 
single particle of an ideal classical gas, and as a result, 
this relationship is true only for a classical gas. 

For a quantum (Bose or Fermi) gas, above probability 
can be given (considering a grand canonical ensemble) 
by Pi = fii/N , where h, is average number of particles 
occupying the ith eigenstate and N{= JA hf) is the total 
average number of particles. Average number of parti¬ 
cles, occupying the ith eigenstate, is, of course, given by 
the Bose-Einstein (—) or Fermi-Dirac (+) statistics 

Ui = e Ci-n)/kT ~ x C 1 ) 

where p is the chemical potential of the quantum (Bose 
or Fermi) gas. Difference between the probabilities (pi) 
of the classical and the quantum gas stems essentially 
from the nonzero fugacity (z = e M / fcT ) of the quantum 
gas. 

Bose-Einstein or Fermi-Dirac statistics as expressed 
above is applicable only in the thermodynamic limit 
where microcanonical, canonical and grand canonical en¬ 
sembles of statistical mechanics reproduce the same re¬ 
sult [8, 18]. However, grand canonical ensemble is conve- 
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nient to work with, in particular, for the quantum (Bose 
and Fermi) gases. 

In the following, we will see how a nonzero fugacity 
plays a vital role in the deviation of Ae 2 from kT 2 c£ l \ 
in particular, for 3-D free [19] and harmonically trapped 
[20] ideal Bose and Fermi gases. 

II. ENERGY FLUCTUATIONS 


temperature dependence of the chemical potentials were 
done by Biswas et al from the above two equations [21], 
They obtained temperature dependent approximate for¬ 
mulas of p not only for the free Bose and Fermi gases, but 
also for the harmonically trapped Bose and Fermi gases. 
Using their appropriate formulas of p, we plot Ae 2 s in 
units of | (kT) 2 in FIG.la for the free Bose gas, and in 
FIG. lb for the free Fermi gas. 


A. For free Bose and Fermi gases 


B. For harmonically trapped Bose and Fermi gases 


For a free gas, we can replace the single particle en¬ 
ergy eigenstates {*} by the single particle momenta {p} 
so that £i —>■ £ P = p 2 /2m where m is the mass of 
the single particle. Degeneracy of this level, within 
the semiclassical approximation, is given by p h yt P , 
where V is the volume of the system. In the thermo¬ 
dynamic limit (A —> oo,V —> oo ,N/V = n = const.), 
straightforward textbook level calculation, for the en¬ 
ergy fluctuation (Ae 2 = £A ef [f^/A] - (]>A G[hi/A]) 2 = 

fo°° 4^P 2( iP-W&ir e P n P p 2 dp] 2 ) of a 3-D 
free Bose gas, results in 


For harmonically trapped case, all the particles are 3- 
D harmonic oscillators, and the energy levels are given 
by £ i = (§+ where oj is the angular frequency of 

oscillations. Although the degeneracy (g t ) of this level 
is i 2 /2 + 3z/2 + 1 [22], yet in the thermodynamic limit 
(A —^ oo, oj —y 0 & Aw 3 = const.), only the first term of 
the degeneracy contributes significantly. In this limit, the 
zero point energy can also be neglected. Now, straight¬ 
forward textbook level calculation, for the energy fluctu¬ 
ation (Ae 2 = A J o °° e 2 n(e)g(e)de - ^r[/ 0 °° en(e)g(e)de] 2 ) 
where g(e) = e 2 /2 in the thermodynamic limit of a 3-D 
harmonically trapped Bose gas, results in 


Ae 2 = 


n„rr\1 ri5 Ll 7/2(z) _ 9 Ll 5 / 2 (z) -| 
V ) 14 Li3/2(2) 4 Li 3/2( z ) J 


for tA > 1 


( 2 ) 


Ae 2 = 


( kT ) 2 [ 12 § H £) 3 - 4 §(£) 6 ] fOT £ < 1 


J C(3) 

) Lis(z) 


( fcT )2[ 12 LgM^ 9 Hf(£)l 


for 7j=r > 1 


( 6 ) 


where T c = ^^-{-^f/ 2 )) 2 '' 3 [19] is the Bose-Einstein con¬ 
densation temperature and Li j(z) = z + fy + + ••■ 

is a polylog function of order j. Similar straightforward 
calculation, for a 3-D free Fermi gas, results in 


Ae 2 


(kT) 2 


~15 Li 7/2 (— z) 
. 4 Li 3 / 2 (— 2 ) 


9Li 5 / 2 ( z) 
^ Li 3 / 2 (— 2 ). 


(3) 


Comparing the above results with the textbook results 
for kT 2 times c v [19], one can easily check that Ae 2 ^ 
kT 2 c v except in the classical limit (z —» 0). To visualize 
how the above results deviate from their unique classi¬ 
cal value (| (kT) 2 ), we can plot the right hand sides of 
Eqns.(2) and (3) in units of | (kT) 2 for the entire range 
of temperature. But, plotting of the same is not an easy 
job until one manages to get temperature dependence of 
z or of p from the implicit relations [19] 


/ Li 3/2 (^) \ 2/3 T c 
V C(3/2) ) T 


(4) 


for 3-D free Bose gas, and [19] 

(-Li 3/2 (-^)r(5/2)) 2/3 = ^: (5) 


for 3-D free Fermi gas whose Fermi temperature is given 
by Tp = (67t 2 7t.) 2 / 3 . Approximate analyses of the 


where T c = )^ 3 [20] is the Bose-Einstein con¬ 

densation temperature for the trapped system. Simi¬ 
lar straightforward calculation, for a 3-D harmonically 
trapped Fermi gas, results in 


Ae 2 


(kT) 2 


-1 r, Li 5 ( z) 
Li 3 (-z) 


Li 2 (-z) ' 

Lif(—-). ' 


(7) 


Once again, plotting Ae 2 , for the entire range of temper¬ 
ature, is not an easy job until one manages to get the 
temperature dependence of 2 or of p from the implicit 
relations 


Li 3 (z) \ 1/3 T c 
C(3) ) T 


( 8 ) 


for 3-D harmonically trapped ideal Bose gas [20], and 

(-Li 3 (-4T(4)) 1/3 = ^ (9) 


for 3-D harmonically trapped ideal Fermi gas whose 
Fermi temperature is given by Tp = ^ (r(4) A) ^ [23]. 
Biswas et al obtained approximate temperature depen¬ 
dent formulas of p from the above two equations [21]. 
Using their approximate formulas of p, we plot Ae 2 s in 
units of 3 (kT) 2 in FIG lc for the harmonically trapped 
Bose gas, and in FIG Id for the harmonically trapped 
Fermi gas. 
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(a): For free Bose gas 

& er/kT 2 c„ (d) 



(c): For trapped Bose gas 

A ,?/kT 2 Cll < d > 



(b): For free Fermi gas 

A <?/kT 2 c v (d > 



0 12 3 4 

(d): For trapped Fermi gas 

A e/kTV d > 



0.0 0.5 1.0 1.5 2.0 2.5 


FIG. 1: Solid lines in FIGs la, lb, lc, and Id (corresponding 
to the approximate chemical potentials in Ref.[21]) respec¬ 
tively represent right hand sides of Eqns.(2) & (3) in units of 
|(fcT) 2 , and Eqns. (6) & (7) in units of 3 (kT) 2 . Dotted lines 
represent classical results. Points represent exact graphical 
solutions. 


For 2-D harmonically trapped Bose gas 

A 2/kT 2 c v ld) 



FIG. 2: Solid line represents Ae 2 in units of 2(kT) 2 . Dot¬ 
ted lines represent classical results. Points represent exact 
graphical solutions. 


III. EXISTENCE OF ‘HUMPS’ IN BOSE 
SYSTEMS 

It is interesting to note in FIG. lc that A e 2 /kT 2 c^ of 
the 3-D harmonically trapped ideal Bose gas is smooth 
unlike its c„, and has a hump just below the conden¬ 
sation point. Maximum in the hump physically means 
maximum average deviation of the energy of the system 
from its average energy both in units of average classi¬ 
cal energy. Appearance of such a hump over the classical 
limit may indicate a discontinuity of c v as well as a phase 
transition around the condensation point (T c ). Here, by 
the phrase “over the classical limit” we mean the average 
deviation of energy is more than the classical average en¬ 
ergy divided by square root of the classical specific heat 
in unit of the Boltzmann constant. 

To measure the energy fluctuation, we need to observe 
energy distribution for different temperatures. From this 
observation, both the average of energy and the disper¬ 
sion of energy can be easily obtained. Temperature can 
be measured from the the width of the distribution. Thus 
one can have experimental data for energy fluctuation 
(dispersion of the energy distribution) for different tem¬ 
peratures. Momentum distribution for the harmonically 
trapped (Bose/Fermi) gas has already been observed by 
applying the technique of releasing the harmonic trap 
[12, 24]. Energy for the harmonically trapped gas has 
also been measured by applying this method [6, 9, 10]. 
Although the specific heat is very difficult to be mea¬ 
sured for the 3-D harmonically trapped Bose gas [17], 
energy fluctuation can be easily measured from energy 
distribution data by applying the technique of releasing 
the harmonic trap. 

That, appearance of a hump in A e 2 /kT 2 c[ cV> over its 


classical limit may indicate a discontinuity of c v , is our 
conjecture. 

Bose-Einstein condensation is possible for 3-D free 
Bose gas. But, there is no discontinuity in its specific 
heat [20]. It is clear from FIG. 1 (a) that, for 3-D free 
Bose gas, there is no hump in A e 2 /kT 2 c^. There is no 
question of phase transition as well as discontinuity of c v 
for an ideal (or weakly interacting) Fermi gas. Obviously, 
there is no hump in A e 2 /kT 2 c£ V> as clear in FIGs lb and 
Id for the ideal Fermi gases. 

Although Bose-Einstein condensation is possible for 2- 
D harmonically trapped ideal Bose gas, there is no dis¬ 
continuity of its specific heat [20]. It is clear from FIG. 
2 that, for the 2-D harmonically trapped ideal Bose gas, 
although there is a hump in Ae 2 / kT 2 ci° l \ it is not over 
the classical limit. 

For the A transition of 4 He, discontinuity of specific 
heat was observed way back in 1935 [25]. Theoretical 
explanation of such discontinuity deals with phonon dis¬ 
persion which effectively simplifies the interacting Bose 
system to a 3-D harmonically trapped ideal Bose gas [26] . 
Considering this simplification, one can easily check that 
a hump over the classical limit must exist in A e 2 /kT 2 Cy 
for the liquid He-4 below the A point. 

Thus we can verify our conjecture for a certain class of 
Bose and Fermi systems. 

Our conjecture linking the maximum in the hump of 
the scaled energy fluctuation (Ae 2 /kT 2 ci°^) and the pos¬ 
sible discontinuity of c v (in unit of cy ) can be justified 
by elaboration through a set of model specific heats as 
shown in FIG. 3 for other Bose and Fermi systems. The 
model specific heat per particle (c„), as represented by 
the dotted line in FIG. 3 (a), is appropriate for a ho¬ 
mogeneous Bose system [19]. We represent average en- 
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(a): For Bose system with a model c v (b): For Bose system with another model c v 

A ,?/kT 2 c„ (d > A d/kT 2 c v (d) 



FIG. 3: The dotted lines in FIGs 3a, 3b, 3c, and 3d rep¬ 
resent model specific heats in units of their respective clas¬ 
sical values (c v /c^) with respect to scaled temperature 
(T/T c ). Dashed and solid lines respectively represent scaled 

energy per particle (e/ \JkT%c£^) and scaled energy fluctu¬ 
ation (Ae 2 /fcT 2 cl cl) ) corresponding to the respective model 
specific heats. 


ergy per particle (e) and approximate energy fluctuation 
per particle (Ae 2 ) in this figure by the dashed and the 
solid lines respectively. We obtain the average energy by 
integrating the model specific heat with respect to the 
temperature (e = f c v dT ), and the scaled energy fluctu¬ 
ation by simply squaring it (Ae 2 = e 2 — e 2 ss ke 2 /cv) so 
that the approximate energy fluctuation becomes exact 
in the classical limit. FIG. 3 (a) represents a special case 
of Bose system with no discontinuity in specific heat at 
T = T c . We do similar exercise in FIG. 3 (b), (c), (d) 
with different model specific heats with finite disconti¬ 
nuity (e.g. in 3-D harmonically trapped Bose gas) at T c 
[1], with infinite (logarithmic divergence) discontinuity 
(e.g. in the A transition of 4 He) at T c [26], and with a 
different finite discontinuity (e.g. in the superconduct- 
ing/superfluid transition of a strongly interacting Fermi 
system) [16], respectively. The humps of Ae 2 / kT 2 Cv s 
over their respective classical limits and the discontinu¬ 
ities of c„s are simultaneously present only in FIG. 3 
(b) and (c) which are appropriate for Bose systems. Al¬ 
though there is a discontinuity of c v in figure FIG. 3 (d) 
which is appropriate for a strongly interacting Fermi gas, 
there is no hump in A e 2 /kT 2 Cv ■ From these figures we 
can also essentially verify our conjecture, and say that, 
for any thermodynamic system, maximum in the hump 
of A e 2 /kT 2 ci cl ^ over its classical limit corresponds to a 
finite or infinite discontinuity of its specific heat c v . 


IV. CONCLUSION 


Here we have analytically explored quantum to classi¬ 
cal crossover in energy fluctuations mainly for free and 
harmonically trapped quantum gases. We have illus¬ 
trated how specific heats differ from energy fluctuations 
for the entire range of temperature. In this regard, we 
have conjectured about possibility of discontinuity of spe¬ 
cific heat in terms of scaled energy fluctuation. 

It is impossible to get exact temperature dependent 
formulas of chemical potentials for the four systems (3-D 
free & harmonically trapped ideal Bose and Fermi gases), 
as because inverses of polylog functions in Eqn.(4), (5), 
(8) and (9) do not exist in closed forms. For this reason, 
we have used approximate temperature dependent for¬ 
mulas of chemical potentials from Ref. [21] . Consequently, 
plots of Ae 2 s in FIG 1, become approximate except in 
the region 0 < T < T c for the Bose gases. To show 
how good our approximate results are, we have com¬ 
pared them with their exact graphical solutions in the 
same figures. In the same figures, we also have plotted 
the classical results just to show the deviations of Ae 2 s 
from kT 2 ci c ^ s. These deviations clearly illustrate quan¬ 
tum to classical crossover in energy fluctuation (Ae 2 ) of 
the four systems of our main interest. 

Although the approximate values of Ae 2 s, for the Bose 
and Fermi systems, match well with their exact graphical 
solutions yet it should be mentioned, that, for the Fermi 
systems, plotting of Ae 2 s had difficulties, specially in the 
low temperature regime, for rapidly oscillatory nature 
of polylog functions in them. For this reason, we have 
plotted the Sommerfeld’s asymptotic forms in FIGs lb 
and Id only for T/Tp < 0.2 [19, 27]. 

It is interesting to note that the scaled energy fluc¬ 
tuation (Ae 2 /fcT 2 c[ cl ' ) ) of the 3-D harmonically trapped 
Bose gas has a hump just below T c . Appearance of such a 
hump over the classical limit may indicate a discontinuity 
of c v as well as a phase transition around the condensa¬ 
tion point. That, appearance of a hump in A t 2 /kT 2 c^ 
over its classical limit might point to a discontinuity of c v , 
is our conjecture for any thermodynamic system. We al¬ 
ready verified this conjecture for a number Bose systems. 
For ideal Fermi gases, as shown in FIG. 1 (b) and 1 (d), 
existence of the hump and discontinuity of specific heat 
is impossible. On the contrary, such possibilities exist 
for a number of interacting or noninteracting Bose gases. 
Although the specific heat is very difficult to be observed 
for the 3-D harmonically trapped Bose gas [17], energy 
fluctuation can be easily measured apparently by apply¬ 
ing the technique of releasing the harmonic trap [12, 28]. 
Thus one can predict the discontinuity of specific heat 
from the existence of the hump over the classical limit. 

By c v , we not only mean specific heat at constant (ef¬ 
fective [29] or exact) volume (V) but also at constant (av¬ 
erage or exact) number of particles N [30]. It should be 
mentioned, that, volume of a harmonically trapped sys- 
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tem cannot be precisely defined although it can be effec¬ 
tively defined for ultracold situation as V e s ~ ( h/muj ) 3 / 2 
[29]. For this reason, c v , for trapped system, is very often 
denoted by C/y. 

In reality, harmonic traps are not isotropic. In that 
case, the angular frequency (w) in all our results should 
be replaced by the geometric mean of uj x , u> y and w,. We 
have not considered spin degeneracy of particles. Consid¬ 
eration of the same would not change any of our results 
as because they are intensive variables. 

The relation between energy fluctuation and specific 
heat is a generic property of microcanonical, canonical 
and grand canonical ensembles of statistical mechanics 
only in the thermodynamic limit which has been taken 
into consideration for the entire calculation of this paper 
[18] . That is why we have expressed both the energy fluc¬ 
tuation and the specific heat in per-particle form. For a 
finite system, this relationship must be different for dif¬ 
ferent ensembles [18]. Extension of our result for finite 
system of different statistical ensembles would be an in¬ 
teresting problem. 

To conclude, this paper illustrates quantum to classical 
crossover in energy fluctuations for a number of Bose and 
Fermi systems. Appearance of the hump in the scaled 
energy fluctuation (A e 2 /kT 2 ) over its classical limit 
may indicate a discontinuity of c v , is our conjecture for 
any (Bose or Fermi) system. Inverse of this statement 


may not always be true as depicted in Fig. 3 (d). All 
the calculations in this article have been done keeping in 
mind the scope of general science readers. Extension of 
our calculations for other spatial dimensions and inter¬ 
acting cases is pretty straightforward and can be taken 
up by the interested readers. They can verify our conjec¬ 
ture for those cases. Proving our conjecture of course, is 
an entirely nontrivial and open problem. 


Acknowledgment 

S. Biswas acknowledges financial support form the De¬ 
partment of Science and Technology [DST], Govt, of 
India under the INSPIRE Faculty Award Scheme [No. 
IFA-13 PH-70]. For an initial part of this work, S. 
Biswas further acknowledges hospitality of the Depart¬ 
ment of Physics, University of Calcutta and financial 
support form the University Grants Commission [UGC] 
under the DSKPDF Scheme. We also acknowledge use¬ 
ful discussions with J.K. Bhattacharjee, HRI, India dur¬ 
ing FIP-2014 conference at the University of Hyder¬ 
abad. We thank the reviewer for his/her thorough review 
and highly appreciate all the comments and suggestions, 
which significantly contributed to improving the quality 
of the paper. 


[1] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999) 

[2[ I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008) 

[3] S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. 
Phys. 80, 1215 (2008) 

[4] A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalat- 
tore, Rev. Mod. Phys. 83, 863 (2011) 

[5] K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van 
Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, 
Phys. Rev. Lett. 75, 3969 (1995) 

[6] J. R. Ensher, D. S. Jin, M. R. Matthews, C. E. Wieman, 
and E. A. Cornell, Phys. Rev. Lett. 77, 4984 (1996) 

[7] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and T. 
Pfau, Phys. Rev. Lett. 94, 160401 (2005) 

[8] S. Biswas, Phys. Lett. A 372, 1574 (2008) 

[9] B. DeMarco, D. S. Jin, Science 285, 1703 (1999) 

[10] J. Kinast, A. Turlapov, J.E. Thomas, Q. Chen, J. Stajic, 
K. Levin, Science 307, 1296 (2005) 

[11] S. Biswas, R. K. Manna, and D. Jana, Eur. Phys. J. D 
66, 217 (2012) 

[12] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. 
Wieman, and E. A. Cornell, Science 269, 198 (1995) 

[13] A. G. Truscott, K. E. Strecker, W. I. McAlexander, G. 
B. Partridge, R. G. Hulet, Science 291, 2570 (2001) 

[14] L. Luo, J. E. Thomas, J. Low Tem. Phys. 154, 1 (2009) 

[15] S. Biswas, Eur. Phys. J. D 55, 653 (2009) 

[16] M. J. H. Ku, A. T. Sommer, L. W. Cheuk, and M. W. 


Zwierlein, Science 335, 563 (2012) 

[17] R. F. Shiozaki, G. D. Telles, P. Castilho, F. J. Poveda- 
Cuevas, S. R. Muniz, G. Roati, V. Romero-Rochin, and 
V. S. Bagnato, Phys. Rev. A 90, 043640 (2014) 

[18] P. Navez, D. Bitouk, M. Gajda, Z. Idziaszek, and K. 
Rzazewski, Phys. Rev. Lett. 79, 1789 (1997) 

[19] R. K. Pathria, Statistical Mechanics, Butterworth- 
Heinemann, 2nd ed. (1996) 

[20] L. P. Pitaevskii and S. Stringari, Bose-Einstein Conden¬ 
sation, Oxford Univ. Press (2003) 

[21] S. Biswas and D. Jana, Eur. J. Phys. 33, 1527 (2012) 

[22] J. L. Powell and B. Crasemann, Quantum Mechanics, 
Addison-Wesley (1961) 

[23] D. A. Butts and D. S. Rokhsar, Phys. Rev. A 55, 4346 
(1997) 

[24] C. A. Regal, M. Greiner, S. Giorgini, M. Holland, and 
D. S. .Jin, Phys. Rev. Lett. 95, 250404 (2005); Q. Chen, 
C. A. Regal, D. S. Jin, and K. Levin, Phys. Rev. A 74, 
011601(R) (2006) 

[25] W. H. Keesom and A. P. Keesom, Physica 2, 557 (1935) 

[26] R. P. Feynman, Statistical Mechanics, Westview Press, 
Colorado (1972) 

[27] Sommerfeld’s asymptotic form of Ae 2 s in units of (fcTp) 2 
are AL + ^(T/Tf) 2 + ... for the free Fermi gas, and 
Aj + /j-(T/Tf) 2 + ... for the harmonically trapped Fermi 
gas. See Ref. [19] for Sommerfeld expansion. 

[28] W. Ketterle, How are Temperatures Close to Absolute 
Zero Achieved and Measured? (Scientific American, 



6 


2004) 

[29] S. Biswas, Pramana - J. Phys. 79, 319 (2012) 

[30] Here “average number of particles” is applicable for 
grand canonical ensemble, and “exact number of parti¬ 


cles” is applicable for microcanonical and canonical en¬ 
sembles. Specific heat per particle is unique for all the 
ensembles in the thermodynamic limit. 


